Ca 2 Y 2 Cu 5 Oio: the first frustrated quasi-lD ferromagnet close to criticality 
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Ca2Y2Cu50io is build up from edge-shared CuCU plaquettes forming spin chains. From inelastic 
neutron scattering data we extract an in-chain nearest neighbor exchange Ji ~ — 170 K and the 
frustrating next neighbor Ji ~ 32 K interactions, both significantly larger than previous estimates. 
The ratio a = | J2/J1 1 « 0.19 places the system very close to the critical point ot c = 0.25 of the J1-J2 
chain, but in the ferromagnetic regime. We establish that the vicinity to criticality only marginally 
affects the dispersion and coherence of the elementary spin-wave-like magnetic excitations, but 
instead results in a dramatic T-dependence of high-energy Zhang-Rice singlet excitation intensities. 

PACS numbers: 75.30.Ds, 78.70.Nx, 74.72. Jt 



Frustrated low-dimensional magnets serve as breeding 
grounds for novel and exotic quantum many-body effects. 
Ca 2 Y 2 Cu50io (CYCO) and the closely related Li 2 Cu0 2 
(LCO) are considered candidates for this type of uncon- 
ventional and challenging physics [l], 0] • These systems 
are in the family of frustrated edge-shared chain cuprates 
(ESC) and their magnetic excitation spectra, as probed 
by inelastic neutron scattering (INS), indicate striking 
puzzles. It was claimed that the dispersion of the mag- 
netic excitations in CYCO shows an anomalous double 
branch 0, Q while LCO exhibits a single, but weakly 
dispersing modeQ. These observations would point at a 
strong deviation of the dispersion from standard linear 
spin wave theory (LSWT), in any realistic ESC parame- 
ter regime. This motivates scenarios with more sophisti- 
cated many-body physics, e.g. due to the presence of anti- 
ferromagnetic (AFM) interchain couplings (IC), causing 
the branch doubling in CYCO 0, However, such a 
scenario invoking strong quantum effects seems to be at 
odds with the observed large and almost saturated mag- 
netic moments ~ 0.9/x B at T < T N = 29.5 K [3, 13], which 
suggests strongly suppressed quantum fluctuations. To 
resolve the situation it is essential to identify the precise 
values of the exchange interactions in these ESC's, both 
within and between the spin chains. To this end it is key 
to measure and at the same time calculate the elementary 
magnetic excitations, ideally for directions of momentum 
transfer in which the excitations depend most sensitively 
on the strength of the in-chain couplings. From scat- 
tering along the a-axis of CYCO, which does not fulfill 
the latter condition, a moderate value of the ferromag- 
netic (FM) nearest neighbor (NN) coupling J\ ps —93 K 
has been extracted with a tiny frustrating AFM next- 
nearest neighbor (NNN) exchange J 2 pa 4.7 K (see Fig. 



1). From a theoretical point of view this is rather unex- 
pected for the ESC chain geometry due to the presence 
of sizable 0-0 2p hopping along the chains. A recent 
reassessment of the exchange strengths based on INS on 
isotopically clean 7 Li 2 Cu0 2 [f| has revealed a relatively 
large FM coupling |Ji| > 200 K, which is more than a 
factor 2 larger than earlier theoretical values 0|. Given 
the structural similarities to CYCO, one expects a larger 
value Ji here too. In fact, high-T 89 Y-NMR data on 
CYCO appear difficult to reconcile with small | J\ | 's [f| . 
Here we show that indeed by measuring with INS the 
magnetic excitations in CYCO along Q = (H, 0,1.5), a 
direction where they are little affected by inter-chain cou- 
plings, one extracts a substantial in-chain J\ pa — 170 K 
and a frustrating J 2 pa 32 K, so that a = |J 2 /Ji| ~ 0.19. 
This indicates an exceptional position of CYCO within 
the ESC family: close to the critical point (CP) of the 
J1-J2 model a c = 1/4, but on the FM side of its phase di- 
agram and in contrast with Li 2 ZrCu04 (ev= 0.3 @) and 
LCO (a = 0.33 @) which are on the spiral side of the 



Ca 2 Y 2 Cu5O 10 (CYCO) 




Figure 1: (Color) Schematic view of the structure of the CuC>2 
chain layer and the main exchange paths of CYCO (see text). 



2 



critical point. We compare the obtained J's to a realis- 
tic 5-band extended Hubbard pd model and L(S)DA+t/ 
calculations, which are in good agreement. The result- 
ing magnetic excitations calculated with exact diagonal- 
ization compare well to the ones obtained with LSWT, 
implying that the coherence of the elementary, spin-wave 
like, magnetic excitations are marginally affected by a 
and quantum fluctuations. However, the relatively large 
Ji and a values obtained affect the thermodynamics [lCj . 
The vicinity to the critical point is probed by e = a — a c 
strongly affects both the magnitude and de- or increasing 
T-dependence of the Zhang-Rice singlet (ZRS) excitation 
intensity for e > and e < 0, respectively. This is mani- 
fest in resonant inelastic x-ray scattering (RIXS), EELS, 
and optics [nj measurements as we will show. 

CYCO has edge-shared CuC>2 chains along the a-axis. 
The Cu 2+ spins are aligned FM along the a-axis (see 
Fig. 1). ac-planes with CuC>2 chains alternate along the 
6-axis with magnetically inactive cationic planes contain- 
ing incommensurate and partially disordered CaY-chains 
which produce a non-ideal geometry in the Cu02 chains. 
These mutual structural peculiarities might be responsi- 
ble for the puzzling strong damping at large transferred 
neutron momenta [H, [i| to be addressed elsewhere. 

Our INS-study has been performed with a fixed final 
neutron energy of 14.7 meV on a 3-axis neutron spec- 
trometer TAS-2 installed at the JRR-3 by the Japan 
Atomic Energy Agency. To analyze the dispersion of 
the magnetic excitations we adopt the model given in 
Ref. [J Then, CYCO has the following main couplings 
J(R), R = (xa, yb, zc): NN and NNN couplings along the 
chain J(1,0,0) = Ji, J(2,0,0) = J2, and the interchain 
coupling (IC) J(0.5, 0,0.5) = J' ic , J(1.5, 0,0.5) = J lc , 
.7(0,1,0) = J b =-0.06 meV, and .7(0.5, 0.5,0) = J ab =- 
0.03 meV. For the small interplane FM couplings Jf, and 
J a b we adopt the values from Ref. [l]. Their contribution 
to the inchain dispersion is negligible. Within the LSWT 
the magnetic excitations dispersion is given by Eq. (2) of 
RefJ: w 2 (q) = A\-B\, A q = J q - J + J„- D, S q = 
J q , where ,7 q = (1/2) ^ r J r exp (zq • r) is the Fourier 
transform of intrasublattice interactions and analogously 
for the intersublattice interactions J q . The dispersion 
along (0, 0, L) shown in Fig. 3 (c) of Ref. H depends only 
on J s = J' ic + Ji C . Its value, as well as the anisotropy pa- 
rameter D, may be found from the INS data for q = 
(0,0,0), (0,0, 1.5): J 2 = ±[ W 2 (0,0,1.5)- W 2 (0,0,0)], 
D = 2 J s - lj(0, 0, 1.5). Using w(0, 0, 1.5) = 5.03 ± 0.03, 
and w(0, 0, 0) = 1.63±0.01 meV we obtain J s w 2.35 meV, 
and D rj —0.27 meV [II] which is very close to J s = 
2.24 meV, and D = -0.26 meV found in Ref. % The 
dispersion along the line (77,0, 1.5) depends only on the 
inchain J's (if the tiny interplane J a b = —0.03 meV is 
ignored). It reads w(q) = Aq = J q — Jo + u(0, 0,1.5). 
These J's may be accessed from the dispersion along this 
line with much higher precision than from the previously 
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Figure 2: (Color online) Dispersion along three lines of the 
first Brillouin zone parallel to the a-axis, obtained from con- 
stant q scans. The LSWT-fit was refined only for the disper- 
sion along (H,0,1.5) and it is shown by thin lines. 



reported data along (77, 0, 0) and we have [12 1 



Ji = -14.69 ± 0.5 (4%) meV w — 170.4 K, (1) 
J 2 = 2.78 ± 0.2(7.6%) meV w 32.2 K. (2) 

The dispersions along (77, 0, 0) (reported in Ref. Q]) and 
(77, 0, 1.25) (given here) are affected by J' ic and Ji C . As we 
have determined only their sum, we adopt J[ c j Jj c = t« 
4/9 for the sake of concreteness (suggested by our band 
structure calculations). Notice that our r- value differs 
from 2 adopted in Ref. [l|. It might be refined empirically, 
if one measures also along (H=1/6,K,L) for any K value. 

With the aim to detect quantum effects beyond the 
LSWT, we calculated the dynamical structure factor 
S(u>, q) using exact diagonalizations (see Fig. 3). In Fig. [5] 
the INS data together with the refined new LSWT-fit are 
shown. S(uj, q) for our set and that of Ref. 1 are shown 
in Fig. 3. The peak positions always nicely follow the 
LSWT-curves, however our set gives a better description 
of the INS-data than that in Ref. [l] (Fig. 4 therein) where 
the artificial double branching was ascribed to AFM IC. 
Indeed, it induces some intensity apart from the LSWT- 
curve, but these intensities are far too weak to be con- 
sidered as a branch doubling. Notice the inflection point 
at 7r/2 for finite a. The total dispersion width is given 
solely by 2| Ji|. 

As an application of our data, we consider the ID- 
magnetic susceptibility x(T) in the isotropic limit (see 
Fig. 4). Despite the uncertainty caused by the unknown 
background susceptibility Xo> the relatively large value of 
| Ji| doesn't allow to extract directly &cw from a l/x(T) 
plot using only data up to 300 K. Instead a much broader 
T-interval up to about 800 K would be required to reach 
the asymptotic high-T limit necessary for a proper quasi- 
linear behavior. Alternatively, higher orders in the high- 
T expansion can be applied [l3[. Hence, the reported 
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AFM values Q C w 15 K0] or weak FM ones Qcw ^ 

+8 K [l5j] are rather artificial. Using the J's from our 
INS-fits we predict instead a markedly larger FM value 



G 



cw 



= 0.5(| Ji | - J 2 ) - J' lc - J ic - J ab = +43.4 K. (3) 



Since we found a < a Cl we readily predict the value 
of the 2D saturation field H s which is here determined 
solely by the total AFM IC like for LCO 0: gH s = 
A(J' tc + J lc ) = -77 A T w 64.8 T, refining an estimate of 
70 T for H || b and g = 2.39 from low- field magnetization 
data |l5j ). Next we consider the magnetic moment in the 
ordered state at low T. Within the LSWT the reduction 
due-to quantum fluctuations is about 6.8% which yields 
1.07/ie to be compared with the experimental value of 
0.92±0.08^b 01 which however is affected by the chem- 
ical reduction effect since about 0.22±0.04 of the local 
moment resides on the O 2p orbitals. 

The exchange coupling strengths can also be deter- 
mined by DFT+C/ calculations. For this we used the 
full potential scheme FPLO [17j (vers. fplo9.01) and 



performed super cell calculations for different collinear 
spin arrangements applying LDA and GGA function- 
al [3, [l^|. The Coulomb repulsion U$d was varied in 
the physical relevant range from 5 to 8eV for a fixed 
J?,d = 1 eV. In our calculations the incommensurate crys- 
tal structure of CYCO can be treated only approxima- 
tively. Thus, we neglect (i) the modulation of the Cu-0 
distances within the Cu02 chains and its buckling and 
(ii) the incommensurability of the Cu02 and the CaY 
subsystems. In particular, the Cu02 chains were treated 
as ideal planar chains reflecting an averaged Cu-0 dis- 
tance of 1.92 A and a Cu-O-Cu bond angle of 94.5°. Fur- 
thermore, we modelled the CaY layer by a Na layer to 
preserve the half filling of the system. The structure of 
the simplified model systems is given in Ref. [5(1 These 
structural simplifications allow a reliable modelling of J\ 
yielding an FM value: ~ -150 K [2(j. In previous studies 
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Figure 3: (Color) Magnetic dynamical structure factor S(ui, q) 
for the J-set of Ref. [2 (left) and for our set (right) from exact 
diagonalizations with L — 14 x 2 and 15 x 2 adopting r = 
D — and Ji C = 2.24 meV, (see Ref. 0. Red line: dispersion 
from LSWT for the two parameter sets (see Fig. 2 and text). 
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Figure 4: (Color) Spin susceptibility \(T) within the isotropic 
ID J1-J2 model using the transfer matrix renormalization 
group method, xo - background susceptibility. 



of closely related ESC [2jJ, |22j for the effects of (i) chain 
buckling and (ii) the cation related crystal fields it was 
shown that J\ is rather robust whereas J2 is strongly re- 
duced by a factor 2 to 3. Thus, our — Ji ~ 150 K is 
considered as a rather reliable lower estimate by about 
10 to 20% with respect to the buckled chain geometry 
in CYCO. However, in view of the drastic dependence 
of J r 2 on these parameters (21, 22 1, a derivation of a re- 
liable value from the applied model structure is difficult. 
Thus, CYCO fits the general experience of a sizable FM 
Ji-value for ESC in contrast to the assignments of only 
a few K, proposed for LiVCu0 4 [1| and NaCu 2 2 [Hf. 
Such small J\s would put them in a region of strong 
quantum fluctuations, harboring the difficulty that the 
observed pitch angle cannot be described classically [ID] . 

Whereas the vicinity of CYCO to the quantum critical 
point only weakly affects the dispersion and coherence 
of the elementary, spin-wave like, magnetic excitations, 
we will show that the amplitude to excite Zhang-Rice 
singlets (ZRS) at typical high-energies, probed by spec- 
troscopic means depends strongly on the frustrating J2. 
and temperature. We have also performed exact diago- 
nalization calculations using an extended 5-band Cu 3d 
02p Hubbard model for CYCO with a standard param- 
eters (Hi. To fit the INS-derived value of J 2 =32 K, the 
tpx,px = 0.59 eV has been slightly reduced as compared 
with LCO (0.84 eV, J 2 = 76 to 66 K) which simulates 
probably the deviations from the ideal chain geometry. 
Mapping the spin states of the 5-band Hubbard model 
onto a frustrated spin model, we obtain J\ = —177.5 K 
and J2 = 32.3 K in full accord with the LSWT- analysis 
of the INS data given above. We stress that the value of 
Ji is mainly determined by the direct FM exchange cou- 
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Figure 5: (Color) T-dependent O-K RlXS-spectrum for xx- 
polarization from a CU4O10 cluster within exact diagonaliza- 
tion using a Lorentzian broadening (half width F = 0.13 eV). 



pling K p d = 65 meV and not by the Hund's rule coupling 



on the O sites J v 



0-5 {U p - U PxPy ] 



0.6 eV as often 



adopted. The significant value of J\ is generic for ESC 
with a Cu-O-Cu bond angle ip < 96° at variance from 
the case of CuGe03 with ip rj 98° causing an AFM J\ . 

As a modern spectroscopy, RIXS provides valuable in- 
sights into the correlated orbital and electronic structure 
(for a review see Ref. 1261 ). Therefore, we also studied 
the T-dependent O-K edge RlXS-spectra for a CU4O10 
cluster within exact diagonalization (see Fig. 5). We find 
a strong decrease of the intensity for the ZRS excitons 
with decreasing T, which is qualitatively in accord with 
general considerations for EELS and optics [llj]. The not 
yet assigned feature observed for CYCO at 300 K near 
527 eV (Fig. 5 in Ref. H3, counted from the 530 eV ex- 
citation energy) corresponds just to these excitations. It 
perfectly agrees with 2.95 eV obtained here. Notice that 



when one sets t t 



0, thereby strongly suppressing 



J2 oc t px px , no ZRS exciton is observed even at 300 K. 
Hence, the frustrated FM differs qualitatively from a pure 
FM in its high-energy response. 

To summarize, we have shown that Ca 2 Y2Cu 5 Oio is a 
frustrated quasi-lD ferromagnet close to criticality. This 
edge-sharing spin-chain compound has pronounced FM 
correlations in the presence of a sizable in-chain frus- 
tration. The main intensity of magnetic excitations in 
S(uj,q) is reasonably well described within LSWT. The 
signatures of the sizable in-chain frustration found here 
cause (i) a characteristic curvature of the in-chain dis- 
persion of magnetic excitations and (ii) Zhang-Rice sin- 
glet features at ~ 3 eV that are strongly T-dependent 
and are detectable by spectroscopies. The spin align- 
ment between the chains below Tv is supported by a 
specific AFM inter-chain coupling which determines the 
saturation field. It causes interesting quantum effects 
that go beyond a linear spin wave description. These 



effects, and the interplay of J' ic and J{ C , provide an inter- 
esting problem not yet investigated in full detail. How- 
ever, such quantum effects are neither strong enough to 
cause a breakdown of LSWT, nor to induce an additional 
branch of magnetic excitations, as suggested previously. 
Our results can further aid in the correct assignment of 
the frustration [23, 24 1 in other ESC, including multifer- 
roics, but especially for the complex chain-ladder system 
(La,Sr,Ca)i4Cu2404i (2^| which also possibly contains 
frustrated FM Cu02 chains as suggested by the edge- 
sharing geometry and sizable NNN transfer integrals |30j 
giving rise to significant AFM J2's. Its nonideal chains, 
as with CYCO, challenge one to look for more sophis- 
ticated but yet solvable theoretical models that include 
incommensurate and disorder effects. 
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In the present EPAPS supplementary online material we present (i) exact diagonalization studies of the dynamical structure 
factor for coupled chains with various interchain couplings of the type as realized in CYCO. (ii) We show our prediction and 
details of the calculated optical conductivity u(uj) (iii) We demonstrate the influence of interchain coupling and spin anisotropy 
on the magnetic susceptibility, (iv) we provide the reader with further details of the LSDA+H7 and GGA+H7 calculations, (v) 
we illustrate the influence of inchain frustration on the dispersion of magnetic excitations focussing on its curvature behavior 
and in particular on the positions of inflection points within linear spin-wave theory. 



The spin-Hamiltonian with uniaxial anisotropy and CYCO geometry 



In general, spin waves are the quanta of small oscillations of spins around the classical ground state of the spin- 
Hamiltonian. So, for the same Hamiltonian, the form of the dispersion may be radically different for different values 
of parameters, since the classical ground state may be different from the quantum one. The spin-Hamiltonian for 
CYCO reads 



H — Ha + Hb + Hab 
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where m enumerates the sites in one sublattice, (r = ±na, n = 1, 2 . . .) determines the neighboring sites within the 
chain, a being the lattice vector along the chain, vector f connects the sites of different chains in the ac plane, R 
connets the sites in different planes, but in the same sublattice in CYCO. We have allowed for an uniaxial anisotropy 
of the exchange interactions. 



(i 



Linear spin-wave theory 

For the collinear AFM the classical ground state is the Neel state, the spins on the A sublattice are directed up, 
and down in the B sublattice. We introduce two different sets of spin-deviation operators 
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The Neel state (AFM ordering of FM chains in the ac plane) is the vacuum state for the operators b |Neel) = 0. So, 
the operator b(a) annihilates the spin deviation from the Neel order (which means a spin-flip for s = 1/2) on the 
sublattice B(A). The Hamiltonian (|S1[) can be rewritten as 



H 
H 



Ho + Hi n t ! 



(S10) 



Bq (a q &_ q + a q 6L q 



A a = S 



^ Jr exp (iqr) - ^ J? + Y Jf 
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(Sll) 
(S12) 
(S13) 



The Fourier transform of e.g. 6 q reads 6 q = y / 2/]V^ mgB exp(— zqm)& m ; N denotes the total number of sites. The 
transverse part of H (|S1|) defines the one- particle hoppings in H (JSTTJ), the Ising part contributes on-site energy 
values. The terms, which contain more than two spin- wave operators enter Hi„ t . The magnons a q , /3 q in an AFM are 
introduced as a mean- field solution of Eq. (|S10[) [? ], which neglects H int . 
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(S15) 
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cosh(9 q &_ q + sinh 6> q a q 



tanh26> q = Bq/Aq. 



Now we specify the expressions for Aq, _B q , which result from the geometry of CYCO and the spin value s = 1/2. In 
the summations over r,R, f we retain only the following terms. For the in-chain exchanges we retain J±, J2, which 
correspond to r — a, 2a, respectively. For the interchain interactions we retain J' ic , Ji C , Jf,, J a b, which correspond to 
f = (a + c) /2, f = (3a + c) /2, R b = b, R a6 = (a + b) /2. Then we find 



Aq = J\ (cos q a a — 1) + J2, (cos2<7 a a — 1) 
+ Jb (cos qbb - 1) + 



q a a q b b 
2 Jab I cos — cos — - 1 



Bq = 2 



J, v cos — — h J ic COS — — 



2J S - D 

(f 



cos 



(S17) 
(S18) 
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Note that interplane interactions enter Aq as they represent intra-sublattice interactions. The value of the gap at 
q = is 



A = ^-Sg, (S19) 
the anisotropy enters the dispersion via the integral value D, which is defined as 

A - Bo = D = ]T ( J f 2 -Jf)-J2 ^ - Jr ) ' ( S2 °) 

f r,R 

and is related to the gap value by the relations 



A = y/D (D + 2B ), 



D = JB* + A^ - B . (S21) 



Dispersion of magnetic excitations and inchain frustration 

Ignoring tiny interactions like Jb for qc = 2(m + 1)tt, m = 0, 1, 2, ... the dispersion along (H, 0,1-5) within the LSWT 
reads (see Eq. 2 in Ref. [H): 

1 — cos x + a (cos 2x — 1) + 
2/3 + 7 (cos3x - 1) + ., 



Ji 



(S22) 

where x — qa, 7 = J3/ | J\ | denotes the third neighbor inchain coupling. The dimensionless interchain coupling 
constant slightly "renormalized" by the anisotropy parameter D: /3d = (J s — 0.5D)/ | Ji |sa 0.17 has been omitted 
in Fig. SI since it doesn't affect the curvature of the dispersion curve in the case under consideration. According 
to our results of the mapping of the 5-band extended Hubbard pd-model on a frustrated spin model 7 is usually FM 
but very small ~ 10~ 2 to 10~ 3 . In many cases it can be therefore neglected. The whole width of the dispersion W is 
given by the odd Fourier coefficients, only: W = 2(| J\ \ +J3 + ...). In the limit i<1 Eq. (S22) can be rewritten as 

1- 4a -97 2 l-16a-8±7 4 6 



Wx - 20* ^ + Otf) . (S23) 

Thus, the inchain frustration a > reduces the quadratic term (which vanishes approaching the critical point, see 
Fig. SI) and strongly affects the quartic term which is negative for the unfrustrated FM (a = 0) but changes it sign 
for a sizable value of a > (1 — 8l7)/16 ~ 1/16 to 1/8, well realized for CYCO with a = 0.19. Thus, a positive sign of 
the quartic term gives direct evidence for a pronounced frustrated ferromagnet. Just the opposite behavior is realized 
near the maximum of the dispersion near the BZ boundary. There the expansion in terms of e = tt — qa reads 

U 2^2(1-7)- 1 + 4 ^ 9 V + 

+ 1 + 1 %- 81 V + 0(s»), (S24) 



Ji 



the inchain frustration a enhances the negative quadratic term (see also Fig. SI) In the case J3 = approximately 
fulfilled for almost all ELC, the inflection point qi v of the dispersion curve (Eq. (S22)) is given by 

-1 ( l- VI + 128a 2 ^ 
q ip a = cos I — 1 . (S25) 



It provides also a clear measure of the inchain frustration a = J2/ | J\ \ present in the system (see Fig. SI). To 
summarize, analyzing the full dispersion curve along the line (H,0,1.5) the presence of the inchain frustration a is 
easily detected. 
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However, we admit that in many cases where the full dispersion curve is not available, a fit of a significant part of 
the total curve might provide a more precise value of a or J2. In such cases the extraction of a from the T-dependence 
of the Zhang-Rice exciton intensity shown in the main text might provide a reasonable alternative. 

Now we briefly consider the (H,0,0) scattering direction where the interchain coupling (IC) is involved. Compared 
with the former "ID "-case we show in Fig. S2 that the IC has a sizable influence on the shape and on the curvature 
of the dispersion. Note that a small quadratic term exists only for finite anisotropy — D > 0. 



A 2 + 4a) /3d + 2/300 + ^)]x 2 + 0(z 4 ), (S26) 



\Ji\ V u 2' 

where 5 = —D/\J\\, (3 — J s j \ J\ |, and /?2 = J%c/\ I J\ I- In the isotropic limit D or 5 — 0, i.e. when the spin gap 
Ao = yj 45/3 + 5 2 vanishes, the dispersion becomes linear (compare Fig. 3 (left) in the main text). 

Like in the former case one observes that the in-chain frustration a strongly enhances the higher order terms near 
the r point. But the interchain coupling reduces this enhancement. Under such circumstances the extraction of the 
IC and a from a a scattering direction affected by both types of interactions is unconvenient. 



Magnetic excitations within the model of coupled chains 

We fix the sum of the two interchain couplings J s = J' ic + Jj C = 2.241 meV at the value derived from the INS data 
and change only their ratio r. First we switch off the frustrating NNN AFM inchain coupling J2 and adopt for J\ 
the value suggested in Ref. [l| (see Fig. 5) One observes roughly two kinds of dispersing peaks: one with a dominant 
intensity well discribed by the LSWT and second one which much less intensity. The latter type of curves has been 
interpreted as "evidence" for an additional anomalous second ME branch in Refs. [l|. A similar behavior is observed 
for all cases with a dominant NN interchain coupling J' ic . In the opposite case of a dominant Ji C that second curve 
almost invisible. For our parameter set that minority peaks are much less pronounced. Thus we conclude, that this 
quantum effect depends on the details of the two main AFM interchain couplings and becomes weaker with increasing 
I J\ I in contrast to what has been suggested previously. Anyhow, a systematic study including the examination of 
finite size scaling is left for a future investigation since it is of less relevance for CYCO. 



The reduction of the magnetic moment at T = 

In contrast to the case of FM ordering, the vacuum state for the antiferro-magnons a q ,/3 q (|S16[) does not coincide 
with the classical Neel ground state |Neel) 7^ |0). Thus, the vacuum |0) contains an finite number of spin deviations 
(|S5|) . (|S7f even at T = 0, (0| h m |0) 7^ 0. Expressing the spin deviation operators a q via a q ,/3 q (|S16p . we obtain 

(0|aL«m|0) = ^^(0|a qflq |0) 



q q v q 



(S27) 



Note, that in our geometry, the summation over q runs over the magnetic Brillouin zone (BZ) which in the present 

b 



case coincides with the lattice BZ — x ^ x — . Finally, for the sublattice magnetization at T = 0, we have 



The optical conductivity within the 5-band extended Hubbard model 



Focusing on the Zhang-Rice exciton, we note that a similar picture as in the RIXS spectra is shown in the predicted 
optical conductivity and EELS spectra (not shown here) . Some insight in the corresponding transitions can be gained 
from the zoomed figure with an artificial small broadening to resolve these transitions (see Fig. S5). 
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site 


x/a 


y/b 


z/c 


Cu 


3/8 





1/4 


On 


7/8 





1/4 


Cu 


1/8 


1/2 


1/4 


Cu 


5/8 


1/2 


1/4 





3/8 





.6270 





7/8 





.6270 





1/8 


1/2 


.6270 





5/8 


1/2 


.6270 





3/8 





.8730 





7/8 





.8730 





1/8 


1/2 


.8730 





5/8 


1/2 


.8730 


Na 


1/8 


1/4 





Na 


5/8 


1/4 





Na 


1/8 


1/4 


1/2 


Na 


5/8 


1/4 


1/2 



Table SI: Crystal structure of the commensurate approximate effective Na2Cu02 compound for DFT+C7 calculation (enlarged 
unit cell and reduced symmetry to allow different spin configuration). Space group P2/M (SG 10), a — 5.6306 A, b = 6.286 A 
and c = 10.5775 A. 



Additional information on the L(S)DA+f7 calculations 

To model the mutually incommensurate crystal structure of CYCO we neglect (i) the modulation of the Cu-0 
distances within the CuC>2 chains and (ii) the incommensurability of the CuO*2 and CaY subsystems as a good 
approximate to estimate the NN exchange along the CuC>2 chains. The crystal structure data for this simplified 
model structure are given in Tab. S2. To calculate the NN exchange Ji for the approximated crystal structure we 
constructed a super cell by the doubling of the unit cell and the reduction of its symmetry to allow for different spin 
configurations (see Tab. SI). The obtained FM NN exchange J\ depending on the Coulomb repulsion U%d and the 
specific functional are depicted in Fig. IS6I J\ depends only weakly on these parameters. The given NN exchange 
J!=-150K is the average between the LDA and GGA result at [73^=6.5 eV. 



Aspects of the magnetic susceptibility 

In addition to ID susceptibilities, for the isotropic as well as the easy-axis anisotropic case calculated using the 
transfer matrix renormaliztion group theory method, also the 3D case has been examined treating the adopted isotropic 
interchain coupling (IC) within the RPA (random phase approximation) (see Fig. S7 and Eq. (S29)). 

X 3D (T) * - *: D(T) m • (S29) 
1 + kxm(T) 



site 


x/a 


y/b 


z/c 


Cu 











O 








0.623 


Na 


1/2 


1/4 


1/4 



Table S2: Crystal structure of the commensurate approximate effective Na2Cu02 structure as a starting model for CYCO. 
Space group FMMM (SG 69), a = 2.8153 A, b = 6.286 A and c = 10.5775 A. 
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Figure SI: (Color) The influence of the frustration a on the dispersion of magnetic excitation for a direction not affected by 
the interchain coupling, e.g. (H,0,1.5) (see Fig. 2 of the main text). Uper: position of the inflection point. Blue curve Eq. 
(S25). lower: the full dispersion for various a- values. 
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2.5 




wave vector qa/rc 



Figure S2: (Color) The influence of the frustration a on the dispersion of magnetic excitation for a direction affected by the 
interchain coupling (IC) e.g. (H,0,0) (see Figs. 2 and 3 of the main text). Black curves: no inchain frustration, red curves: 
inchain frustration included. Broken curves curves: Ji C included, but we adopt J; c =0 for the sake of simplicity. 
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Figure S3: The dynamical structure factor from exact diagonalizations of two coupled chains with 14 sites for each chain at 
various energies u and momenta q. Red line: projected maximum position of S(ui,q) very close to the LSWT result. Blue 
dashed line: intensity from the weak minority peaks and the inchain parameters proposed in Ref. 1. 




co (meV) 



Figure S4: (Color) The same as in Fig. S3 for the new parameter set given in the main text. 
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Energy (eV) 

Figure S5: (Color) The optical conductivity <j(uj) for a CugOi2 cluster with periodic boundary conditions within the 5-band 
extended Hubbard Cu 3d 02p model and parameters given in the main text. Upper: Broad energy region and spectrum 
broadened by V = 0.3 eV. Notice the NN Zhang-Rice exciton visible in the spectrum at 300 K and its lacking at T — 
(compare Fig. 5. in the main text). Lower: The same as in the upper for the region of Zhang- Rice singlet transitions. The 
spectrum has been broadened by T = 0.05 eV, only, to make all transitions visible. Notice the main Zhang-Rice singlet exciton 
near 3 eV. 
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Figure S6: (Color online) Calculated FM NN exchange Ji as function of the Coulomb repulsion U^d- 
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Figure S7: (Color) Left: spin susceptibility x{T) fitted within the ID J1-J2 model supplemented with isotropic AFM IC in 
2D treated within the RPA . The latter is measured by the parameter k — 2 (J' ic + Ji C ) / \ Ji \. Right: The same as left for an 
anisotropic easy-axis inchain coupling J%. 



